A Remark on Integrability of Stochastic Systems 
Solvable by Matrix Product Ansatz 

V. Karimipour [] 



Department of Physics, Sharif University of Technology, 
P.O.Box 11365-9161, Tehran, Iran 
Institute for Studies in Theoretical Physics and Mathematics, 
P.O.Box 19395-5746, Tehran, Iran 



Abstract 

Within the Matrix Product Formalism we have already introduced a multi-species 
exclusion process in which different particles hop with different rates and fast 

particles stochastically overtake slow ones. In this letter we show that on an open 
chain, the master equation of this process can be exactly solved via the coordinate 
Bethe ansatz. It is shown that the N-body S-matrix of this process is factorized into 
a product of two-body S-matrices, which in turn satisfy the quantum Yang-Baxter 
equation (QYBE). This solution is to our knowledge, a new solution of QYBE. 
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1 Introduction 



The static Matrix Product Ansatz (MPA) as formulated in ||, relates the problem of finding 
the steady state of a one dimensional homogenous stochastic system , to a quadratic algebra 
and its representations. 

Looking back at the history of the subject, one can look at static MPA and its final formu- 
lation by Krebs and Sandow||, as an abstractization of some recursion relations between 
certain quantities for systems of different sizes. In this way the technique of using recursion 
relations [f§]fl acquires a purely algebraic character, and the results obtained by the for- 
mer method are rederived || by computing matrix elements of product of operators. 
Conceptually MPA is the generalization of the ordinary Bernoulli measure where matrix 
elements of operators replace numbers to produce correlations. 

Many of the questions regarding the steady state and the correlation functions can then be 
answered at least in principle, by purely algebraic manipulations. We say "in principle", 
since in general, finding the representations of the algebra may be as difficult as the original 
problem. The special cases where the static MPA has been used for calculating concrete 
physical quantities include the one-species asymmetric simple exclusion process (ASEP) on 
a ring and on an open chain with and without an impurity ||,@„[f§,|J,[]lO|](see also |Tl|, [p^| 
and references therein). The works on the multi-species asep include [O [pE|,Pfl, and M. 



When one comes to the dynamical MPA [15] | 16fl , the corresponding algebra becomes time- 
dependent and much more involved. Even when all the particles are identical, only in certain 
limited cases the dynamical algebra has been used for calculating physical quantities, namely 
the symmetric exclusion process [plj and a non-equilibrium model of spin relaxation [|16|]. In 
these two cases the algebra is simple enough to be used for finding the spectrum of the cor- 
responding quantum Hamiltonian and r-point density correlators. These two cases happen 
to be integrable. In fact they belong to a 10-parameter family of diffusion-reaction systems 
which are equivalent through a similarity transformation to a generalized Heisenberg chain 



T6| , which is known to be integrable. 



For a general stochastic process, dynamical MPA, like the static one is just a rewriting of 
the master equation and its real power is revealed only when one can solve the resulting 
algebraic relations or can find a representation of the algebra. 

It is useful to agree upon a notion, however vague, of MPA solvability. Let us call a process 
static-MPA (dynamic-MPA) solvable, if a nontrivial representation of the static (dynamic) 
MPA algebra can be found, or the ensuing recursion relations can be solved. 
At present there is no definite connection between static or dynamic MPA-solvability and 
integrability, although from the above examples it appears that the dynamic-MPA solvable 
systems are integrable, since for these models the entire spectrum has been found. However 
it is not known what kind of integrable systems are MPA solvable. 
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To see if there is any connection between MPA-solvability and integrability, it is desirable to 
study a number of examples which hint to such a connection and to investigate the common 
properties of these examples. 

In this letter we provide an example of a multi-species stochastic process which on the one 
hand is related to a very simple algebra in the MPA formalism and on the other hand is 
integrable in the sense that its N-particle S-matrix is factorized into a product of two-body 
S-matrices, which in turn satisfy the quantum Yang-Baxter equation (QYBE). The algebra 
is generated by p + 1 generators E and D^i— 1 • • -p with relations JTJ 

DiE = -Di + E i = l...p (1) 

= : j d 1 -v 1 d 1 

Vi - Vj 

where Vi < v% < • • • < v p are finite non-zero real numbers. In this process a particle of type 
j hops with rate v j and when it encounters a particle of type i, with Vi < Vj the two particles 
interchange their sites, with a rate Vj — Vi as if the fast particle stochastically overtakes the 
slow one. This model seems to be natural as a simple model of one-way traffic flow. We 
have shown elsewhere Q,0 that on a finite chain this process is MPA-solvable in the above 
mentioned sense. We now consider the same process on an open chain and show that the 
master equation of this process can be solved via the coordiante Bethe ansatz. That is, by 
a suitable transformation we turn the master equation into an eigenvalue equation, and find 
the eigenvalues and eigenfunctions of this equation, the latter being in the form of Bethe 
wavefunctions. In this case, since we have p species of particles, the two particle S-matrix 
is not a c- number as in [17], [IS] but is a p 2 x p 2 matrix. Thus the factorizability of the 
N-particle S-matrix puts a stringent condition on this two-body S-matrix, in the form of the 
Yang-Baxter equation. We will show that this is indeed the case, and as a byproduct provide 
a new solution of QYBE with spectral parameter. 



2 The Master Equation 

Suppose there are N particles of different types on an infinite one dimensional chain. The 
basic objects that we are interested in, are denoted by P ai ,a 2 ,-a N ( x ii x 2, • • - %n) which are 
the probabilities of finding at time t, a particle of type a\ at site X\, a particle of type a 2 at 



site X2, etc. Following [TJ[], these functions to define probabilities only in the physical region 
X\ < x 2 < ■ ■ ■ < Xj\f. The hypersurfaces where any of the two adjacent coordinates are 
equal are the boundaries of the physical region. In the subset of the physical region where 
Xi + i — Xi > 1, Vi we have free hopping of particles and the master equation is written as 
d 



2 



— ("qi 4" ' ' ' fa JV )-Pa 1 .02,---ajv( a 'l) ^2) " " ' > %N'i t) • (3) 

We assume hereafter for convenience that the time variable has been rescaled so that all the 
hopping rates are dimensionless quantities, i.e. they are ratios of actual hopping rates to 
some standard hopping rate. 

Extending the above equation to the whole physical region, where Xi + % — Xj > 1, causes 
terms on the boundary surfaces to appear on the right hand side of the master eqution, 
which should be fixed by a choice for the boundary condition. The choice of boundary 
condition in fact determines the type of interaction between particles (see ||19|| ). On any 



such surface where for example Xj + i — x, = 1 for some i, we supplement the above equation 
with the following boundary condition, where we suppress for simplicity the time variable 
and all the other coordinates: 

vpP a ,p{x, x) = v a P a ^(x, x + 1) + (vp - v a )P^ a (x, x + 1) (3 > a (4) 

Pp, a {x, x) = P /3)Q (x, x + 1) (3 > a (5) 

The master equation (0) and the boundary conditions (§,|5]), replace the very large number of 
equations which one should write by considering the multitude of cases which arise according 
to which group of particles are adjacent to each other. Instead of giving a general proof we 
provide a few examples in the two and three particle sectors. 
Consider the two particle sector. From the definition of the process we have: 

d 

—Pi 2 (x, X + 1) = ViP 12 (x - 1, x + 1) + (v 2 - Vi)P 2 l(x, x + 1) - V 2 Pl2(x, x + 1) (6) 
ot 

d 

— P 2 i(x, x + 1) = v 2 P 21 (x - 1, x + 1) - (vi + {v 2 - wi))P 2 i(z, x + 1) (7) 

These are exactly the equations which are obtained from combination of and (^, |5|). 
As a couple of examples in the three particle sector, consider P 123 (x,x + l,x + 2) and 
P 2 i3(x, x + 1, x + 2). From the definition of the process we have 

d 

Pi 2 3{x,x + l,x + 2) = viPi 23 (x - l,x + l,x + 2) + (v 2 - vi)P 2 i 3 (x,x + l,x + 2) 



dt 



+ (v 3 - v 2 )P 132 (x, x + 1, x + 2) - v 3 Pi 23 (x, x + 1, x + 2) 



d 

P 2 i 3 (x, x + 1, x + 2) = v 2 P 2 i 3 (x - l,x + l,x + 2) + (v 3 - vi)P 23 i(x,x + l,x + 2) 



Ot 



(v 3 + (v 2 -vi))P 213 (x,x + l,x + 2) (9) 



These are precisely obtained from combination of (|3|) and |5|). This pattern repeats in all 
sectors. 
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3 The Bethe Ansatz Solution 

We rewrite P au0 t2,-a N ( x i, %%■,"•■> x n\ t) as 

P ai ,-a N {x u ...,x N -t) = e^ E+v ^+- v ^ t (v ai ) xl ■ ■ ■ {v aN ) x ^ ai ,.. aN {x u ■■■,x N ;t) (10) 
and insert it into (0) to obtain 



(Xx — 1, X2, • • • , Xn', t) + ^ ai ,a 2 ,-a N (, x li x 2 ~ 1; ' ' ' ; x N'i t) 



+ ■ ■ ■^a 1 ,a 2 ,-a N {xi, x 2, ■ ■ ■ ,X N - l;t) = -E^ ai ,a 2 ,-a N ( x l , x 2 , ■ ■ ■ ,X N ;t) (11) 

In terms of \P the boundary conditions ([|, [5]) are rewritten as 

v 

^ a> p{x,x) = v a ^ a ^{x,x + 1) + — (vp - v a )Vp ta {x,x + 1) (3 > a (12) 



^3,a(x, X) = V a ^ Jx, X + l) 



f3>a 



(13) 



Hereafter we use a compact notation as follows. ^ is a tensor of rank N, whose components 



are \I/ 



Cri,cr2,-"CW 



^i, x 2 , • • • , xjv)- The boundary conditions are written in the form 

*(■••, e,e,---) = b M+ i*(---,e,e + i,---) 



(14) 



where b&^+i is the embedding of b (the matrix derived from (JT^, |I3"D ) in the locations A; and 
fc + 1. 

b fc ,fc+i = l®---®b®---®l (15) 
In the two species case b is equal to 



/ 



Vi 



. \ 



V2 J 



In the p-species case it is : 



y ■ 



(16) 



(17) 



where (Eij) k> i = S ik 5ji and % = iW(ij)- 

To solve the eigenvalue equation (|TT]) , we write \I/ as a Bethe wave function 

*(x 1 ,---,x w )=^A ae ^- x . (18) 

Here x and p denote n-tuples of coordinates and momenta respectively, and the summation 
runs over all the elements a of the permutation group. For each element a of the permutation 
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group, the corresponding coefficient (a tensor of rank N) is denoted by A a . Inserting fll8|) in 
(El), yields 

A a e ia{pyx (Yl e ~ l<T[P3) + E) = (19) 

a 3 

From this, one obtains the eigenvalues (as one can remove a from the summations in the 
parenthesis) 

E = -(e~ ipi + e~ ip2 + ■■■ e~ ipN ). (20) 

The Bethe wavefunction should also satisfy the boundary condition fll^p. Inserting QT8D in 
([14] ) we obtain 



(21) 



We now note that Xj and £ are arbitrary and the coefficient of £ is symmetric with respect 
to the interchange of p k and Pk+i- This interchange is effected by the element a k of the 
permutation group. Thus symmetrizing with respect to this interchange, we obtain 



(l - ^^'b M+1 k + - e^b fe)fe+1 W & = 



(22) 



where aa k is the product of the elements a and a k in the permutation group and A aak is the 
corresponding coefficient. Thus we obtain 

A aak = S k>k+1 (a(p k ), a(p k+1 ))A a (23) 

where 

Sfc, fc+1 (zi, z-j) = 1 ® • • • ® S(z u z 2 ) <8> ■ ■ ■ ® 1 (24) 

and 

S ( ^ 1 ,^) = -(l-2 1 b)- 1 (l-2 2 b) (25) 

Here we have denoted e iPk by S(zi, z 2 ) is the two particle S-matrix. The above equation 
allows all the A a J s to be calculated recursively in terms of A\. The first few members of the 
permutation group are 1, <7i, a 2 , (Ti(T 2 , 0201 and CiO^Ci = o"20"io"2- Accordingly we find from 

An = Su{pi,P2)A A a2 = S 23 (p 2 ,p 3 )A 1 A aia2 = S 23 (pi,p 3 )A ai (26) 

Ann. = Sn(pi,P3)A a2 Ann<n = Si 2 (p 2 , p 3 )A aia2 . (27) 

Moreover we have 

Ar 2CTl( T2 = S23 (PhP2) A a2ai (28) 

At this point one may think that the master equation ([3]) can be solved by Bethe ansatz for 
any type of matrix b and hence for any type of nearest neighbour interaction. The answer is 
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that we are not yet done with the Bethe ansatz solution, since for consistency of the solution 
we should check that different ways of obtaining a coefficient yeild the same result. It is 
sufficient to check this consistensy only for two cubic elements of the permutation group 
fl20l| , namely o\<JiG\ and 020\<Ji which are equal as elements of this group. Thus we should 
have 



A„ 



(29) 



This yields via 



(28) the following Yang-Baxter equation: 



Sl2(P2,P3)S23(pi,P3)S 12 (pi,P2) = S 23 (pi , P2) Sl2 (Pi , Pz) S 2 3 (P2 , Ps) 



(30) 



which is crutial for consistency of the method. 

Remark: With the transformation R := PS, where P is the permutation operator, we 
obtain the more familiar form: 



#23(P2,P3)#13(Pl,P3)#12(Pl,P2) = #12 (Pi, P2)#13(Pl, P3)#23(P2, Pa) 



(31) 



We can calculate the S-matrix from ([Tj|) and (25). For example in the two-species case we 
find 



S(z, w) 



v-i z— 1 



1-VlW V! ( \ w-z 

1 — V\W 
VlZ—1 



\ 



1—V2U) 

V2Z 



(32) 



I J 



Note that the last diagonal element is different from the others. We have checked that this 
matrix satisfies exactly equation (|30|) , although its proof requires lenghty calculations which 
we will not reproduce here. It has moreover the following properties 

1- In the special case where V\ = v 2 we have S(z,w) = i=si which is nothing but the 
S-matrix found in [17]] for the single species case. Note that in the latter case the S-matrix 
is a c-number and (|30| ) is satisfied trivially. 

2- S-\z,w) = S{w,z) 

3- S(z, z) — 1 V vi,V2 

One can easily obtain from ( |T8"D and (^), the S-matrix for higher values of p. The final 
result is 



S(z, w) 



P 1-VijW 

/ j ~ 1 ■ 



1,3 V 



En + E 



Vi 



Kj 3 



w — z 



ViZ^ 



Ein <8> E. 



,2 » 



J' 



(33) 



where Vij 



In conclusion we have provided an example of a multi-species stochastic process (the number 
of species being arbitrary), which is both static-MPA solvable (see the introduction) and 
integrable on an infinite chain. Integrability of the model is ensured by testing a very 
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stringent condition on the S-matrix, namely the 2-particle S-matrix satisfies the Quantum 
Yang-Baxter equation. This solution is to our knowledge a new solution of QYBE. Given 
these facts, it would be quite interesting to see if this process is also dynamic-MPA solvable. 
At a concrete level and for sectors of low number of particles, one can use the Bethe wave 
functions to calculate the same quantities that were calculated in |18|], namely the diffusion 
and drift constant, and see what is the effect of overtaking on these characteristics. Questions 
specific to the multispecies case can also be answered. For example, given an initial situation 
of particle 2, I steps behind particle 1, what is the probability that at a given time particle 
2 overtakes particle 1, or what is the average passing time. 
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